Centre de Physique Theorique - CNRS - Luminy, Case 907 

F- 13288 Marseille Cedex 9 - France 

Unite Propre de Recherche 7061 

^ EXACT LYAPUNOV EXPONENT 

a^ : FOR INFINITE PRODUCTS 

OF RANDOM MATRICES 



(N 






R. LIMA and M. RAHIBE|] 



^^ Abstract 

O I In this work, we give a rigorous explicit formula for the Lyapunov expo- 

nent for some binary infinite products of random 2x2 real matrices. All 

"^ ■ these products are constructed using only two types of matrices, A and B, 

J>^' which are chosen according to a stochastic process. The matrix A is singular, 

Y ■ namely its determinant is zero. This formula is derived by using a particular 

O i decomposition for the matrix B, which allows us to write the Lyapunov ex- 

r-| I ponent as a sum of convergent series. Finally, we show with an example that 

O ■ the Lyapunov exponent is a discontinuous function of the given parameter. 
> 
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1) Introduction 

The product of random matrices appears in the study of disordered sys- 
tems [1] as well as in the context of dynamical systems [2]. The Lyapunov 
exponents are one of the tools to study these products. They are related 
to physical quantities in disordered systems [3]. For example, in the thight 
binding model or Anderson model [3], the localisation length of the wave 
function is proportional to the inverse of the Lyapunov exponent. In dy- 
namical systems theory, products of random matrices often arise as a non 
trivial approximation wich mimics strongly chaotic behavior in deterministic 
systems. 

In spite of the important and numerous results obtained in the theory 
of random matrices [4,5,6], there is no general method for calculating the 
Lyapunov exponents [7,8,9], although for a few examples there is an explicit 
formula [10,11]. 

In this paper, we present some examples of products of random matri- 
ces, where we were able to determine the Lyapunov exponents as a sum of 
explicitly convergent series. In all these examples, we deal with infinite bi- 
nary random products, built with only two types of 2 x 2 real matrices, A 
and B which are chosen according to a stochastic Bernoulli or Markovian 
process and one of the matrices, say A, is singular, i.e. its determinant is 
zero. Markovian processes are chosen to mimic the correlations existing in 
dynamical systems exhibiting weak chaos [12]. The key of our study is the 
use of a particular decomposition when the matrix B is non-singular (but if 
B is singular, the calculation can be done directly). The Bernoulli case was 
studied by Pincus [13], who gave assymptotic results for the Lyapunov expo- 
nent, but which do not enable an actual computation, see remark 2.3 below 
(also Derrida and Hilhorst [14] have obtained a similar formula for a partic- 
ular product). Instead, we give an explicit form of the terms of that series. 
In some cases, this permits either an explicit summation of the series, or a 
control of the convergence, therefore leading to approximative results. Fur- 
thermore, we could show that the convergence of the series is exponentially 
fast. 



The Markovian products are treated in the same way as the Bernoulh 



one. 



The paper is organized as follows. In section 2, we derive the general 
formula for the largest Lyapunov exponent 7. Section 3 is devoted to the 
simple case where 5 is a singular matrix. In section 4 we perform the decom- 
position of B wich, in turn, leads to the decompostion of -B" and in section 5 
we derive the expression of 7. In the last section, we analyze the continuity 
of 7 as function of parameters. 

2) Lyapunov exponent, general formula 

We consider A and B, 2x2 real matrices, A and B, where A is a singular 
matrix. Now let Vn : 



V, 



N 



XnXn~1 -^2-^1 



be a binary product of A^ matrices, where the matrix Xj is either A ot B, 
the choice being made by a stochastic process. In the Bernoulli case, we 
have Xi = A with probability p {0 < p < 1) and Xi = B with probability 
q = p — 1, {i > 1). In the Markovian case, the transition probabilities are 
given by 

fPr(X„+i = A/X„ = 5)=pi 

Pr{Xn+i = B/Xn = B) = l-pi = qi 

[ Pr{Xn+i = A/Xn = A) = 1-P2 = q2 



and 



Pr{Xi = A) 
Pr{Xi = B) 



Pi 



Vl +P2 
_P2_ 

Pi +P2 



Po 

go 



with < pi < 1 and < p2 < I- 

By definition, the Lyapunov exponent 7 is 



(2.1) 



(2.2) 



7 = lim T7 log II Vn 



(2.3) 



7 is independent of the choice of the norm 



Since the matrix A is singular, it can be written by a change of basis in 
one of the following form : 



A 














A 









^=(n n) (2-4) 

or, 

A=[; „j (2,5) 

If A is of the type (2.5), the Lyapunov exponent is 7 = —00. This is 
straightforward to show, since A'^ = 0. 

Therefore we suppose, without loss of generality, that A has the form 
(2.4).,.e.A=(^ I). 

/6ii(n) 6i2(n)\ 
If we write B^ in the form B^ = , then the result 

Xhiin) h22{n) I 
obtained by Pincus [13], in the Bernoulli case, follows : 

00 
7 = plog I A \+Y.p\l-pT\og I 6n(n) |, (2.6) 

n=l 

Where p^(l— p)" is the probability to obtain the subproduct AB^^A. Even 
if (2.6) was proved only in the Bernoulli case [13], the same argument extends 
to the Markovian case, leading to the following proposition. 

Proposition 2.1: Let { Vn } be an infinite product of random matri- 
ces satisfying the Markovian distribution law (2.1) and (2.2), where A is a 
singular matrix given by (2.4) and B is general. 

Then 

Pi °° 

7 = ■ log I A I + V PoPipagi""^ log I hii{n) \ (2.7) 

Pi + P2 ^1 

where hii{n) will be explicitly computed in section 4. 



Remark 2.2: The proof of this proposition is analogous to that given in 
[13] once we notice that the probabihty to find the product AB^A is poPi^r"^ 
in the Markovian case whereas it is p^(l — p) in the Bernoulh case. 

Remark 2.3: Notice that, in order to give an exphcit value for 7, from 
(2.6) for the Bernouli case or (2.7) for the markovian case, the key point is 
the computation of 6ii(n), a problem which was not addressed in [13]. This 
is the object of the next sections 

According to (2.7), it is possible to study the Bernoulli and the Markovian 
cases in the same manner, by writing 



7 



Pr{A) log I A I +L ^ x"-^ log I 61 



iW 



n=l 



where 



' Pr{A) = p 
< L = p^(l -p) 



(. X 



p 



in the Bernoulli case, and 



Pr{A) 



Pi 



L 



pIp2 



Pi +P2 



(2. 



(2.9) 



(2.10) 



P1+P2 
(X = qi 

in the Markovian case. 

3) The case of two singular matrices 

In the case where B is singular, it is easy to calculate bii{n). Indeed, as 
explained, above we have two different cases, 
either 

^-«-*(o' o)« 



,-WO At 



^ = « lo oj« 

where Q is an invertible matrix and Q^^ is its inverse. 

In the former case 7 = —00, by the same argument as when A is given 
by (2.5). 

In the latter case, if Q is written as : Q = ( ) , then we have, for 

V q'21 q'22 / 

^ > 1, biiin) = Afe"'giig22- And it is then easy to obtain 7; for a Bernoulli 
product, we have 

7 = 7b =plog I A I +(l-p)log I Xb I +p{l-p)\og\ ^TT^ I (3-1) 

and for a Markovian product, we have 

P^ 1 I \ I , ^2 , I \ I , P1P2 1 I ^11 I fn o\ 

1 = 1M = ^ log A + ■ log Af, + ■ log 3.2 

P1+P2 P1+P2 P1+P2 Tr(fi) 

where 611 denotes the first entry of the matrix 5, and Tr{B) is the trace of 
B. 

In the formulae (2.11) and (2.12), we notice a non-linear term, originating 
in the non-commutativity of the matrices A and B. 

4) Normal form and computation of B^ 

When B is non-singular, we introduce a decomposition, which enables us 
to determine 611 (n) for all positive integers n. 

Let i? be a real, non-singular 2x2 matrix, which we write as 

B =1 det5 I n{-ip)Bn{ip) (4.1) 

where det B is the determinant of B, TZ(ip) = i . | is a matrix 

V sm (/9 cos (p J 

rotation in the plane with an angle if and TZ{—(f) is its inverse. The angle (p 
is determined by 

ta„(2^) = ^H^ (4^2) 

O12 + O2I 

5 



where bij are the entries of B. 

We call the matrix B, the normal form of B. We have four different 
types of normal forms B, depending on the sign of the determinant and on 
the eigenvalues of B, which can be real, non-degenerated eigenvalues, real 
degenerated eigenvalues and conjugate complex eigenvalues. 

We now define the quantity : 



where 



and 



P^\/l|l (4,3) 

k = bi2 sin^(v9) - 621 cos^(v9) + -(611 - 622) sm{2ip) 

I2 = hi2 cos^(v9) - 621 sin^(v9) - -(fen - 622) sin(2v9) 
p will be used in the expression of B. 

Remark 4.1: Replacing B hy C = B, accroding to (2.3), 7 is 

det-B 
shifted by an additive constant, log | det 5 |. Therefore, we can assume 

without loss of generality, that | det -B |= 1- 

We now consider four different cases. 
Case \ : B is hyperbolic symplectic 

The eigenvalues of B are real non-degenerated and have the same sign, 
Ai,2 = eexp ±0" with a non zero and e = 1 or — 1 ; then, we get 

(cosh a p sinh a \ 
- sinh a cosh a j \ ■ ) 

Case II : i? is hyperbolic, non-symplectic 

In this case the eigenvalues of B are also real, non-degenerated but they 
have opposite signs, namely Ai = exp a and A2 = — exp —a or — Ai and — A2. 

6 



Although we have two possible forms for B, we can deduce one from the 
other by changing the sign of ip. We therefore retain only one form for B: 

sinh a p cosh a " 
' - cosh a sinh a ' \ ■ ) 



Case III : B is parabolic 

The eigenvalues of B are given by Ai^2 = e, (e = 1 or —1) and we have 
either 

B = e{\ tV (4.6) 



,0 1 , 
or 

«-(-k ?) (-) 

Notice that we can pass from (4.7) to (4.6) by a rotation in the plane 
with an angle — , which is equivalent to a change of the sign of (p. In the 
following, we define the normal form of a parabolic matrix as given by (4.6). 

Case IV : B is elliptic 

The eigenvalues of B are complex conjugate, Ai^2 = exp±ia, in which 

case 

(cos a — p sin a \ 

1 . . (4.8) 

-smcr coso" / ^ ' 

Therefore, if i? is a non-singular matrix, we can write it in the form given 
by (4.1) where B is given by one of the expressions (4.4), (4.5), (4.6) or (4.8). 

We can now easily compute B"- in each of the previous cases. We sum- 
marize the result in the following proposition. 

Proposition 4.2: Let Bhe a, non-singular matrix which a normal form 
B defined by (4.1) 

7 



Then 

5" = 7^(-v9)i3"7^((^). (4.9) 

If B is symplectic hyperbolic, given by (4.4), then 

(cosh na p sinh na \ 
1 • u u • (4.10) 

— smh na cosh ncr j ^ ^ 

If B is non-symplectic hyperbolic, given by (4.5), then either 

(coshno" psinh.n(j\ 
^ ■ h h \ (4-11) 

- smh na cosh na j ^ ' 

if n is odd; or 

(sinh na p cosh na \ 

^ h ■ h \ (4-12) 

— cosh na smh na j ^ ' 

if n is even. 

If B is parabolic, given by (4.6), then 

i3" = e"('J ""^j^^'V (4.13) 

If B is elliptic, given by (4.8), then 

(cosncT — p sin na\ 
1 . (4.14) 

— smna cos na j ^ ' 

5) Lyapunov Exponent : explicit formulae 

We are now ready to perform the analysis of the formula (2.8) in the case 
where the matrix B is non-singular. 

As previously, without loss of generality, we suppose that | det-B |= 1- 

We now give the value of the largest Liapunov exponent in each of the 
four cases treated in the previous section. 



Case I : B is an hyperbolic symplectic matrix 

For n > 1, we have 

^11 (^) = e"[cosh(no") + cosh(a) sinh(?T,o") sin(2(/9)] (5-1) 

where 

1 1 

cosh a = -(p+ -) (5.2) 

2 p 

In (5.1), we can take o" > 0. Indeed, if a is negative, we change the sign 
of ip and use —a instead of a. 

Notice first that if it exists an integer no > 1 such that 6ii(no) = then 

7 = — oo. 

The condition 6ii(no) = is equivalent to 
cosh a ■ 



sin(2(/9) tanh(nocr) 
and therefore we can construct some products for wich 7 = —00. 

On the contrary, if we suppose that for all integers n > 1, 611 (n) 7^ 0, and 
we define the quantities 

5 = 1 + sin(2v9) cosh(a) (5.3) 

and 



1 — sin(2(/9) cosh(a) 

T = 

1 + sin(2v9) cosh(Q;) 
We will use r and 5 in the expression of 7. 

We now distinguish two cases. 

a) (5 = 0: 
A straithforward calculation gives 



(5.4) 



7 = 7b =plog I A I -(1 -p)(T (5.5) 



for the Lyapunov exponent in the Bernoulh case, and 

7 = 7m = \ log A ■ a (5.6) 

Pi + P2 Pi + P2 

for the Markovian case. 

b) (5^0: 
Depending on whether r is zero or not, we have different expressions for 

7- 

i) r = : we easily obtain the largest Lyapunov exponent for a Bernoulli 
product: 

7 = 7b =plog I A I +(1 -p)o-, (5.7) 

and for a Markovian product: 

7 = 7m = \ log I A I H ■ a (5.8) 

Pi + P2 Pi + P2 

ii) r 7^ : Here we obtain 7 as a convergent series given by 



7 = Pr{A) log I A I +(Y^a 

+ 13^ log f + i: E^=i x^-^ log I 1 + re-2"- 



(5.9) 



where L is defined by (2.9) and (2.10). 

Notice that if 7iv is the sum of the N-th first terms in (5.9) of 7, the error 
is 

I 7 - 77V |<| r I exp(-2(A^ + 1)(t)Lx^+^ (5.10) 

and therefore it is exponentially small since < x < 1. 

Case II : B is an hyperbolic no symplectic matrix 

This case is similar to the previous one with only a slight difference con- 
cerning (5.9). By applying the decomposition given by (4.1) we obtain 

^ l + sin(2v9)cosh(a) ^„^ ^ „ 1 - sin(2v9) cosh(a) ^_„^ 
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Again, we can suppose that a > without loss of generahty. And, as 
above, if there exists an no > 1 such that 6ii(2no) = or 6ii(2no — 1) = 0, 
then 7 = — oo. 

Recall that S and r are defined by (5.3) and (5.4). 

If 5 = or if 5 7^ but r = 0, the largest Lyapunov exponent is given 
by the same expressions as in the previous case: (5.5) and (5.6), or (5.7) and 
(5.8). 

Instead, ii 6 ^ and r 7^ 0, we obtain 



7 =Pr(A)log|A|+(^a+^log|§ 
+LE^=ia:""'log| l + (-l)"re 



00 ™n— 1 ]„„ 1 I /^ 1 \n^^,—2na ^ ' 



This serie is convergent, indeed we have 

I7-7JV |<^a;^+^A (5.13) 

where 

A = max{\ log I 1 + (_i)A^+ie-2(^+i)- ||, | log | 1 + (_i)^+2g-2(A^+2V ||)_ 

(5.14) 

Case III : B is a parabolic matrix 

The normal form of B is 



and thus, for n > 1, 



"-ij \ 



h^^{n) = e"[l + -nbsm{2ip)]. (5.15) 



We suppose that 611 7^ for all integers n > 1. The trivial case where 
sin(2(y9) = corresponds to an infinite product of diagonal and triangular 
matrices. In this case, 

7 = Pr{A) log I A I . 

11 



When, instead, sm{2ip) ^ 0, we have, for a Bernoulh product: 

7 = j9 log I A I +p{l — p) log I '"°^^^'^-' I 
V(l-p)E^=i(l-pr-Mog|n 



(5.16) 



(5.17) 



n6sin(2i^) I 

and, for a Markovian product, 

7 = -^ log I A I +^^ log I '-^^ I 

' P1+P2 ° ' ' P1+P2 ° ' 2 I 

+ P1+P2 Z-n=l U Pi J log I n + ^3i^2v,) I- 

These series are convergent since we have 

2 1 A?" + 1 

~ n6sin(y9) l — x (1 — a;)^ 

and < a; < 1. 

Case IV : B is an elliptic matrix 

In this case we also obtain series for 7, but it is very difficult, in general, to 
evaluate the rest of the corresponding partial sums. Indeed, in this situation 
, as expected, the summability of the series in the expression of 7 is related 
to the arithmetic properties of a, when the latter is irrational {mod 27r). On 
the another hand, if a is rational {mod 27r), 7 may be given as a sum of a 
finite number of terms. 

As above, by using the decomposition (4.1), we obtain 

^11 (^) = cos(no') + sinh(a) sin(no-) sm{2ip) (5.19) 

where 

sinh(a) = ^(p--) (5.20) 

with a G H. 

If 611 (no) = for some positive integer no, then 7 = —00. 

Suppose now that 611 (n) 7^ 0, for all positif integers n, and a is rational 
r 
{mod 27r), i.e. a = -27r where r,s E IN* are irreducible. Thus 

s 

7 = plog I A I +- ^ -£ (1 -p)^log I 6n(j) I (5.21) 

12 



for a Bernoulli product, 
and 

7 = ^^log I A I + j}P\ ^E (1 -Pir'log I &ii(j) I 

P1+P2 [Pi + P2){1 - [l - PlY) ~{ 

(5.22) 
where 

2 IT '2lT 

bii (j) = cos( tt) + sinli a sin(2(/9) sin( vr) (5.23) 

s s 

for a Markovian product. 

6) One parameter family of products: an example. 

One may wander if a limit procedure may permit the computation of 7 
in the elliptic case for irrational a from the corresponding formulae (5.21) or 
(5.23), when cr is rational. 

More generally, the continuity of the largest Lyapunov exponent as a 
function of the amplitude of the disorder in disordered systems is known as 
being an important issue. In this section, we illustrate the results of the 
previous sections by giving an example where, indeed, 7 is a discontinuous 
function of a parameter a. 

We consider the family {Poo(tt) = {A ,B{a) , pi , P2), a G R,^*, pi > , 
P2 > 0} of infinite Markovian products of the singular matrix A = 





and a family B{a) of matrices, depending on a parameter a. pi and p2 are 
the transition probabilities defined by (2.1). 

li Pi + P2 = 1, then we recover the case of a Bernoulli product. 
B{a) is defined by 

where (p = —j, a > is a fixed parameter; and p = e" with a G IR^*. 

The condition 611(a) = gives us 

a = an = —Argch{coth.{ncr)) (6.2) 

13 



and 7(an) = — oo. 

Now for all a G IR"^*, we define 

J \ I l\ > P1+P2 ^ P1+P2 &V 2 / /n o\ 



f ^ En=i (1 - PiT-' log(l + ^i^e-2-o)) 



Since the series in (5.3) is sommable, then /(a) is defined and continuous. 
But since, for a ^ an, we have 

7(a) = f{a), (6.4) 

We obtain 

lim 7(a) = /(a„) 7^ 7(a„) (6.5) 

for all n G IN*. Thus 7(a) is a discontinuous function of a, for each a = an, 
n G IN. 
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